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PREFACE 


The Agriculture and Resources Inventory Surveys Through Aerospace Remote 
Sensing is an 8-year program of research, development, evaluation, and appli- 
cation of aerospace remote sensing for agricultural resources, which began in 
fiscal year 1980. This program is a cooperative effort of the National 
Aeronautics and Space Administration; the U.S. Agency for International 
Development; and the U.S. Departments of Agriculture, Commerce, and the 
Interior. 

The work which is the subject of this document was performed within the Earth 
Resources Research Division, Space and Life Sciences Directorate, at the 
Lyndon B. Johnson Space Center, National Aeronautics and Space Administration. 
Under Contract NAS 9-15800, personnel of Lockheed Engineering and Management 
Services Company, Inc., performed the tasks which contributed to the 
completion of this research. 
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1. INTRODUCTION 


Radar scatterometer measurements were taken by the National Aeronautics and 
Space Administration (NASA) C-130 aircraft at 1500 feet above ground level as 
part of the Agricultural Soil Moisture Experiment (ASME) near Colby, Kansas, on 
July 18, 20, 21, and 22 and on August 8, 9, and 11, 1978. These data are a 
measure of the efficiency of surface backscatterl ng as seen from the sensor. 
Four radar scatterometers were mounted on the aircraft: 0.4, 1.6, 4.75, and 
13.3 GHz. Their location on the aircraft Is Illustrated In figure 1. Each 
sensor collected data by looking aft with Incidence angles between 5” and 50*. 
All angles were sensed simultaneously. 

Soil samples from several layers were collected from preselected fields of 
approximately 40 acres on each of the 7 flight days. These san^les were 
weighed, oven dried, and weighed again so that the moisture content of the 
layers could be calculated. They will be used for comparison with the 
scatterometer data. 

The purpose of this report Is to document a method of converting the scatterom- 
eter computer-compatible tape (CCT) data Into disk files, wherein each file 
contains the date, the scatterometer data, and the ground reference position of 
the data within the sampled field. This conversion Is accomplished by execut- 
ing three programs on the National Advanced Systems AS/3000 computer, located 
In the Earth Observations Division Laboratory (EODL) In Building 17 at the NASA 
Lyndon B. Johnson Space Center (JSC). A listing and discussion of each program 
are given In this document. 


2. REQUIRED INPUT DATA 

As the plane flew down the flight line, the scatterometer data were recorded on 
tape In analog form. The tapes were sent to the Sensor Analysis Laboratory 
(SAL) In JSC Building 15. At the SAL, the analog data were subjected to a 
Doppler frequency shift filtering which separated the data by Incidence angle. 
Then, the data were digitized and manipulated so that all Incidence angles 
representing approximately the same target area were grouped together. 
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Figure 1.- Bottom view of the NASA C-130 aircraft 



In the case of the Colby ASME, the output CCT’s from the SAL contain four files 
for each sensor on each run over a flight line. Only files 1 and 4 are 
required by the programs described In this document. The two files containing 
sensor calibration and bandwidth Information (files 2 and 3» respectively) are 
not uued. 

File 1 contains the aircraft's flight parameters. Each record In the file con- 
tains the roll angle, pitch angle, drift angle, altitude, ground speed, and 
true heading angle of the aircraft referenced by the Auxiliary Data Annotation 
System (ADAS) time of measurement. Mew aircraft flight parameters were 
measured approximately every six-tenths of a second. A listing of a sample 
file Is shown In figure 2. 

File 4 on the CCT contains the scatterometer data. Figure 3 Is a sample listing 
of this file. Each record In the file contains the backscatter coefficients. In 
decibels, for approximately the same area as viewed from 10 different Incidence 
angles. They are referenced by the ADAS time when the aircraft passed over or 
closest to the target area. If the drift angle Is nonzero, the areas viewed at 
various angles are unregistered and decorrelated. 

A great deal of aerial photography was collected at the Colby test site. 
Photographs taken at an altitude of 8000 feet were used to construct controlled 
strip mosaics of each flight line. Additional aerial photography was acquired 
as the aircraft collected scatterometer data at altitudes of 1000 and 1500 feet. 
The acquisition time and the frame number of every photograph were recorded on 
the analog data system. This made It possible to determine the aircraft's posi- 
tion at the frame t1».as. The camera positions and frame numbers were plotted on 
transparent overlays by the JSC Cartographic Technology Laboratory. Additional 
overlays were made showing the location of the sampled fields In each flight 
line. All overlays were made at the same scale as the strip mosaics. 
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Figure 3.- A sample listing of the file containing the scattero6eter data. 
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3. PROGRAM EXPLANATION 


3.1 SCAT 

The first proRram used to process the scatterometer data Is a Fortran program 
called SCAT. A listing of the program and Its execute file are given In the 
appendix. This program reads files 1 and 4 after they have been transferred 
from the CCT to a disk. The program uses these data to compute the Inter- 
section of the aircraft's negative z-axis with the ground In a scene-based 
coordinate system. 

Three Inputs from the terminal are requested by the progrmn. The first Input. 
AMISS. Is the northward distance (In feet) of the plane from the southern field 
boundaries. If the flight line runs east-west, or westward from ^e eastern 
field boundaries for north-south flight lines, at the beginning of the flight 
line. The second Input. YUP. Is the crosstrack distance (In feet) that the 
airplane's position Is from the reference field boundary at the end of the 
flight line. Figure 4 Is a diagram of a flight line and shows the distances 
represented by AMISS and YUP. Both distances are measured (In millimeters) 
with the overlay on the strip mosaic. The final Input reouested Is called 
CODE. It Is a three-symbol numeric Identifier for the day. sensor, and polari- 
zation of the data. The codes used for the Colby ASME are shown In table 1. 

Program SCAT creates two output files. One file contains the downtrack and 
crosstrack locations of the airplane at each time that aircraft flight param- 
eters were available. The second output file contains the scatterometer data 
and the location of the aircraft's negative z-axIs Intersection with the ground 
with respect to the beginning of the flight line. The file also contains the 
distance (In feet) of footprint decorrelation caused by aircraft drift. This 
distance represents the difference In crosstrack location of the 5* and 50* 
sensor footprints. 
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Figure 4.- > sai^ile plot Illustrating how and where to measure the variables MISS and YUP 


TABLE 1.- COOES FOR THE COLBY SCATTEROWETER DATA 


Code 

Day. 1978 

Sensor, GHz 

Polarization* 

« 

199 - 7/18 

13.3 

VV 

2 

201 - 7/20 

4.75 

VH 

3 

202 - 7/21 

1.6 

HV 

4 

203 - 7/22 

0.4 

HH 

5 

220 - 8/08 



6 

221 - 8/09 



7 

223 - 8/11 




*VV • vertical transmit and vertical receive, 

VH ■ vertical transmit and horizontal receive, 

HV « horizontal transmit and vertical receive, and 
HH ■ horizontal transmit and horizontal receive. 


3.2 PLOT 

The second program, PLOT, is a Statistical Analysis System (SAS) program. This 
program reads the two files output by progrmn SCAT and produces a printer plot 
of the aircraft’s ground track along with the corresponding 20* scatterometer 
data. The plot is at the same scale as the strip mosaics. The 20* incidence 
angle ms chosen because it was Judged to be most responsive to row direction 
effects. The advantage of this will be explained later in this paper. A por- 
tion of the plot is shown in figure 5, and a program listing is given in the 
appendix. At this point in the analysis, the scatterometer data are referenced 
by the distance (in feet) down track and the distance from the southern or 
eastern field boundaries, depending upon the direction of the flight line. It 
is necessary to know the ground reference position of the data within the 
sampled fields. This is accomplished by using the overlays in conjunction with 
the plot. Both overlays, one with photographic position and ttie other with 
field boundaries, are placed on the plot in the following manner. First, a 
time is found when the aircraft’s flight parameters are available and when an 
aerial photograph was taken. The time represented by each asterisk in the 
flight path plot can be found by using the exclmnation points plotted along 
side. The exclamation points are time marks when the aircraft clock was at the 
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Figure 5.- A portion of the plot depicting the aircraft's ground track 


minute or a multiple of 10 seconds after the minute. Next, the overlays are 
placed on the plot so that the photographic position from the overlay Is on top 
of the asterisk representing the same time. The solid line paralleling the 
plot of the flight path should be even with the southern field boundaries for 
an east-west flight or with the eastern field boundaries for a north-south 
flight. Figure 6 Is an Illustration of the plot with the overlays in place. 

The downtrack distance (In feet) of each field's closest boundary to the begin- 
ning of the flight line Is read from the plot. These distances, along with the 
dimensions of the corresponding fields, are written in a separate file. This 
file Is used as an Input to the next program. 

When the procedure Is carried out In the manner described above, a discrepancy 
may become apparent. If the overlay and plot were registered at an extreme end 
of the flight line, the overlay and plot positions will not necessarily corre- 
spond at the opposite end. This Is due to a lack of sufficient accuracy In 
recording the airplane's flight parameters. Therefore, It Is recommended that 
the overlay and plot be registered at several points along the flight line. 

The plot of the 20* scatterometer values is useful for a final, small 
adjustment in the alongtrack direction. 

The row direction of the fields can be determined from the aerial photography. 
Since the 20® scatterometer angle Is highly responsive to row direction 
effects, knowledge of how the row direction changes from field to field can aid 
in determining exactly which points fall in the samoled field, 

3.3 GRID 

The GRID program reads in the boundary data file, along with the file contain- 
ing the scatterometer data that was created by program SCAT. A listing of this 
Fortran program and Its execute file are given in the appendix. A separate 
output file Is generated for each sampled field within the flight line. Each 
record In the file contains the distance (In feet) from the 30® footprint to the 
northern and the western field boundaries, along with the corresponding scatter- 
ometer data and decorrelation distance. The output files can then be combined 
In the manner which best suits the analysis techniques that will be used. 
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Figure 6.- Illustration of the plot of the aircraft's ground track with 
overlays indicating the photographic position and field boundaries. 



APPENDIX 


PROGRAM LISTINGS AND EXECUTE FILES 



APPENDIX 


PROGRAM LISTINGS AND EXECUTE FILES 

A listing for program SCAT and its execute file are provided in figures A-1 and 
A-2, respectively. The listing for program PLOT is given in figure A-3. 

Figures A-4 and A-5 are respectively the listing and execute file for program 
GRID. 
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* the distance (In FEET) TO THE 
FIELD BOUNOAfliES FOR THE FLIGhI 
« USED TO PUT time MARKS ON TK' 


toty 

I TEN 

THE OUTPUT PARAMETERS TO UNIT 13 ARE: 


SOUTHERN/EASTERN 

iT Parameters 

fHE PLOT 


? MEAR 
CATX 

SCATY 

KTEN 

ISIGZ 


THE FOOTPRINT DECORRELATION DISTANCE. 

THE OOWNTRACK DISTANCE (IN FEET) OF THE 
SCATTEROMETER DATA 
THE DISTANCE (iN FEET) TO THE SOUTHERN/EASTERN 
FIELD 80UNQAH1ES FOR The SCATTEROMETER DATA 
US60 TO PUT THE TIME MARKS ON THE PLOT 


. _ . _ riME _ . _ 

AN ARRAY THAT CONTAINS THE SCATTEROMETER DATA FOR 
ALL ANGLES IN ONE FIELD 


ccccccccccccccccccccccccccccccccccccccccrcccccccccccccccccccccccccccccc 

icAL^foTlIllolfloTY’Islof ll^ENlslOUAC^Ecliio) fO£CORR(510) 

REAL SCATS(510)*Y0FF(S10) «SCATY(S10) *0RIFT(S10) 

REAL SIGZ(IO) »SCATX(S10) fKTEN(SlO) fTRACK(SlO) *SREEO(310) 

REAL SHEAR (510) 


C 

C 

C 


THESE statements ASK FOR AND READ INFORMATION FROM THE TERMINAL 


IN F3.0*F3.0*) 


C 

C 


WRITE(16*161) 

161 FORMAT (♦ INPUT AMISS AND YUP 
READ (15. til) AMISS* YUP 

151 .f6rmaT(F3.0.F3.0) 

WRITE(16.162) 

162 FORMATCTYPE IN THE THREE NUMBER CODE FOR OAY/FREQ/POLAR* ) 
REA0(15.152)IC1»IC2.1C3 

152 format (311) 

YUP« ( (AMISS'YUP) /25. A) *1000.0 

AMISSs(AMISS/2S. A) *1000.0 

CHGSNsl.O 

YOLO«0.O 

HT0T*0.0 

TOTY(1)»0.0 

TOTX(f)«0.0 

THESE statements READ IN AND MANIPULATE THE AIRCRAFT OATA 

Figure A-1.- Coinputer listing for program SCAT. 


SCA002S0 

SCA00300 

scAooalo 

$CAg0320 

SCA00330 

SCA00340 

SCA003S0 

SCA00360 

SCA00370 

SCA00380 

SCA00390 

t CAOOAOO 
CAOOAIO 
SCA00420 
SCA00A30 
SCAOOAAO 
SCA00450 
SCA00460 
SCA00470 
SCA004G0 
SCA00490 
SCAOOSOO 
SCA00510 
5CA00520 
SCA00S30 
SCA00S40 
SCAOOSSO 
SCA00S60 
SCA00S70 
SCAOOSBO 
SCA00590 
SCA00600 
SCA00610 
SCA00620 
SCA00630 

scaoEgsS 

SCA00660 

SCA00670 

SCA00680 

SCA00690 

SCAOU700 

scAOono 

SCA007|0 

SCA00730 

SCA00740 

SCA00750 

SCAU0760 

SCAU0770 

SCA00780 

SCA00790 


A-2 


OF POOR 


qUW.»TY 


FILCi SCAT 


FORTRAM A 


CONVCRSATIONAL HONITOR SYSTEM 


1F(ILINE.E0.1>CH0|N«-I.d 

!F(fLfNE.CO.*)CHQSN«-f.5 
IF(tLINE,C0.5)CHeSNa-I*0 
DO 15 H«l,2 

r£ad(10*105)0ummy 


101 


ENOsZOl ) ET« iHiNtSCC *ROLL'f PITCHfORIFT ( 1 1 «ALTf SPEED ( I > 

12*12«T1S*F 
2»T/6*F6.2) 


RHAT(T2«FS.2*T12*12«TlS*F4.1»T26tF6.2tT36»F6.2«T46«F6.2* 
6tF6. * 


20 


IFLAG»0 
DO 1 I«N510 
REAOdOtlOl* 

F0RMA1 

minSIloatI ImIni 

YOFF ( i ) *ALT*TAmOTOR) 
numbaC»I-1 
lMtN«|MtN«1000 
IseC>lFlX(SEC*10.0) 

Haiplir:!!---'' 

‘:hk.6t,3o)!fla5;i 

- 

^ ]FLA6.EQ.2I ACSEC ( 1 ) »ACSEC ( 1 > «3600.0 
PRACK (n«ORlFT(I)*HCAD 
CONTINUE 
YUP«YUP/(ET*SPEED(NUH6AC) ) 

NNsNUMBAC»l 


iF(IFLA6.E6.inflM£(I)»IT|MF(I)-A00^ 
lF{lFLA6.^0.2nT|M|u)*IT!MElI)*60000 

* TRA( 


THESE STATEMENTS COMPUTE THE AVERAGE FLIGHT DIRECTION 


14 


NN«FLOAT(NN) 

"2t4? 

HAVG 


«HTif*VRACKCK) 
iHTOT/ANN 


THIS SECTION COMPUTES THE LOCATION OF THE AIRCRAFT 

00 11 l«ltNUMBAC 
TOTD»0.0 
DO 18 KJ«i«10 

TOTO*TOTO*ORIFT ( JJ) 

AV6D«( (TOTD/1 0.0) *1.141593) /180.0 
1(1)«(TAN(AV60)*1350.0)*CH6SN 


18 


:ORRI 


, , 

EL?TME«?L0AT < IT IME ( I ♦ 1 ) -I time < 1 > > /lO.O ^ 
ANGLE* ( I (HAVG-TRACK(l) )) *3. 1415931/180.0 


»P»ELPTME*SPE£0(I) 
[SP»OISP*COS(ANGLE) 

I sp«o I sp*s I N 1 angle ) *CHGSN 

^^Y6iFFi{YOFF(?)iYOFF(I-l) )*CH6SN 

T0TX(L)*T0TXCL-1)*X0 




TOTY(L)*TOTY( 
fxiL?*xJlSP 


,)*YU 


liJ’ 

I SR* 


YDIFF* (VUP*TOTX<U-YUP*TOTX lL-1 ) ) 


TO _ 
TOTY(U 
CONTINUE 


YOISP*YOFF(I)*AMISS 


12 
11 

THIS SECTION DETERMINES WHICH AIRCRAFT DATA IS AT 10 SECONDS 
I«2*NUMBAC 

nIi)*q.O . . 

MP«IT!Mf (p/100 

■ iTlMEni-nTEMP^lOO 
FF.LE.3)ITeN7l)*TOT 
■.6T.96)ITE 



8 CONTINUE 


Y(L)*100.0 

EN(I)*TOTY(L)*100.0 


THIS SECTION READS IN THE SIGMA 2EH0*S 


1CA00800 
kAOOOlO 
SCA00020 
SCA00830 
U00840 
[A00850 
[A00860 
^iAOOOTO 
SCA00880 
SCA00890 
SCA00900 
>SCA00910 
SCA00920 
SCA00930 
SCA00940 
SCA00950 


lAOlOOO 
sdAOlOlO 
SCA01020 
SCA0|03C 
SCAO 040 
SCAOIOSO 
SCAO 1060 
SCAO 070 
SCAO 080 
SCA01090 
5CA01100 
SCAOllIo 
SCAO !|o 
SCA01130 
SCA0U4O 

lasr 

SCAO .] 
SCAOI 
SCAO 
SCAO 
SCAO] 
scAoi; 
scAoij 

SCAOll 
SCA012SO 
SCA01260 
SCA01270 
SCAO 
SCAO 
SCAol 
SCAOI 
SCAOI 
SCAOli.. 
SCA01340 
SCA01350 
CAO 1360 
CA01370 
SCA01380 
SCA01390 
SCA01400 
SCA01410 
SCA01420 
SCA01430 
SCA01440 
SCAO 1450 
SCA01460 
SCA01470 
SCA01480 
SCA01490 
SCA01500 
SCA01510 
SCA01520 
SCAO 
SCAO 
SCAO 
SCAO 
SCAO 
SCAO 


Figure A-1.- Continued. 


A-3 


ORIGINAL RAGE IS 
OF POOR QUAtnY. 


riLCt SCAT 


FORTRAN A 


CONVERSATIONAL. HON 1 TOR SYSTEH 


107 


102 



H»l»3 

<li*107)0UNHV 
RMAT(A4) 

LAOaO . 

2 KalfSlO 


17 




DO 17 N«l*10. 

lSl6Z(K*N)aIFlX(S16Z(NI*10.0) 
lUMBSOaF 


NUMBSOaK 
MlNaFLOATiKMIN) 

SCATS (K) a(MlN*60»0) «SEC 
RMl.4aKHIN*1000 
KSECalF!x(SCC*10a0» 
KTIME«K)aKM “ " 
KClKaKTI^F 


IFCKCHK.Gt 
IF(K 


20 


[F(K^LA6. 
FCKFLAG. 
IFCKFLAG, 

CONTI^U^^ 


.30)KFtA6aI 


HK.LT.-S0000)KFLAGa2 
(K)«KTIMJ 


i)KT|ME 

Sid?! 


(K)aKTIHf 

(K)aSCAT! 


IK)-AO 

(K>^60 

(K)*36 


il8;i 


these statements compute which SCATTEROMETER data is at 10 SECONDS 


Ka1,NUMBSO 
;N(X)aC ‘ 


DO 9 , - 

KTEN(i«)a0.0 
KTEMPaKTlME«K>/100 
KOlFFaKTlME(K)-(KTEMP»100» 
MMalSIGZ(K.A) 


! F<KDlFF.LE.3»KTENtK)aFLOAT<MH»/10.0*0.5 
F(KniFF.GT«96)KTEN(K}aFLOATlMH)/10.0«0.S 
F(KTFN(K) •GT.5.0)KTEN(K)aKTENIK)>1.0 
9 CONTINUE 

THESE STATEMENTS FIND OOwNTRaCK LOCATION OF THE SIGMA ZERO*S 




32 


33 


UNUMBSO 
IK)a0.0 
SCATX(K)a0.0 


FC 

[•!?1 


(SCATS<K).LT.ACSeC(p) 
fSCATS(K).LT.ACSEC(i«inGO 


GO 


IfU.GT.NUHBAOGO to 203 
60 TO 32 ... . 

RANGE laACSEC ( I ♦ 1 ) 'ACSEC C 1 1 
PAN6E2aSCATS (K I -ACSlC Tl ) 
RATI0aRANGE2/RAN6El 


TO 3 
TO 33 


Lai 

IF(1.F0.1)60 TO 13 
ALOFFaTOfx fL) -TOtx (L-l) 

YYOFFaTOTY (U -TOTY <L-i) 
OCOFFaOECORR<L»-OeCORRIL-l» . . 

SCATY (K ) a ( YYOFF^HAT 10) ♦TOTY CL-1 ) 

SCATX (K) a(ALOFF»RATlO) ♦TOTX CL-ll 

Smear (k) a ( dcoff*ratio) ♦decorr cl-i ) 

GO TO 3 

SCATX (K) aTOTX (L) *RATIO 
CONTINUE 

THESE statements WRITE OUT THE FILES TO DISK 


203 CONTINUE 

DO 4 IaI«NUMBAC 
4 WRITFtl2*103> 


103 


- 10 . 

1 *NUMBSO 


format (• 

104. FORMATM **31 
1 IX.FA.OtlXt 

t TOP 
NO 


a03>TOTX(I) .TOTYdltlTENm 
•♦3F10.2) 


:i *C2*C3« iLlNE.IRUNf SMEAR (K) .SCATX (K) .SC aTY (K) . 
IGZ(KtML) .MLal.IO) ., . 

i.1X.11.«/'.iI«1X.F4.0.1X»F6.0.1X.FS.O. 

iOU) 


Figure A-1.- Concluded. 


ICA01S70 

scAoIeeo 

SCA01890 

SCA01900 

SCA01910 

SCA0l9|0 

t CA01930 
CA01940 
SCA019S0 
SCA01960 
SCA01970 
SCA019B0 
SCA01990 
SCAU2000 
SCA02010 
SCA02020 
SCA02030 
SCA02040 
SCA020S0 
SCA02060 

18:81818 
SCA02090 
SCA02100 

i§:8j 


SCA02190 


SCA02300 

l8:8Hi8 

8:81318 


A-4 


.H UINAL PAGE IS 
OF POOR QUALITY 


M page is 


PILEt SCAT 


Exec 



if Sill! tl sgifii g !?|5i! 

tl tSatav E IpIrm 

i?iS8 


(PERM HECPM P 




SCAT 


RECPm 

RECPM 


CONVERSATtONAt MONlTOf SYSTEM 



t imniii lii 


Figure A-2.- Execute file for program SCAT. 


A-5 



ORIGINAL PAGE IS 
OF POOR QUAlfTY 

riLCt PLOT SAS A CONVERSATIONAL MONITOR SYSTEM 


mm ?8JK1 81 


then T1C«*I 


OUTPUT PARTI t 
OUTPUT PART2t 


YTIC«i.O«flCI 
■“ X Lg 40000 THEN 

''pUtHjStIV®" 

inpile P1LE4I 

INPUT 1*15 X 6.0 *28 TTIC 4.0 MAS Yl 4.11 

A4«Tiyi*io.o)/3.oT«iooo.oi 

IF TTIC ■ 0.0 THEN fTlC«.l 

Y2>AA*4000.0I 

YVTIC«TT1C*10I 

YTIC2«{YYTIC/3.0) •1000.0*4000.01 
DROP YYTIC AAI 


EN 


OUTPUT PART3t 
OUTPUT PART4f 


40000 TM 
, 20000 TH 
. ’H 

merge parti PART3I 

BY XI 

DATA PLOT2I 

MERGE PART2 PART4I 
8Y XI 

POOC PLOT OATAmPLOTlI 

plot x*Y»**» X*YT!c»»I* X«Y2«»*» 
HAX|S*-3000 TO 4000 BY 1000 
VAXfs«0 TO 40000 BY 1000 
HSPACE»10 VSPACE«fl 
MPPP»0 3000 3150 7333 
VREF»0 to 40000 BY SOOO 
VPOS«350 I 

PROC PLOT OATA«PLOT2l 

PLOT X*Y«*«* X*YTTc»*I* X»Y2«**» 
HAXIS«-3000 to 9000 BY 1000 
VAXfSB20000 TO 62000 BY 1000 


X*YTIC2»» I '/OVERLAY 


X*YTIC2«» I '/OVERLAY 


HSPACE«10 VSPACEbB 
HPPFbO 3000 3150 ‘ 


VP 

VP 




733 


000 TO 62000 BY 5000 
0 I 


Ir 


Figure A-3.- Computer listing for program PLOT. 


A-6 


ORIGINAL PAGE IS 
FILC» 9R10 FORTRAN A OF POOR ^ 


CONVERSATIONAL MONITOR SYSTEM 


^CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC^ 


IS PROGRAM IS the THIRQ OF THREE PROSRAMf USED TO^PROCESS. , 
t COLSY ASmI SCaTTEROmETER DATA, THIS PROOHAM SEPCRaTCS THE 
TA INTO files for EACH SAMPLED FIELD ALONG ThE FLIGHT LINE. 


TH 
THi 
DA 

THE INPUT VARIABLES FROM UNIT 9 AHCt 

JXT FOR explanation 
. T LINE 


e?5e : m jfjisSt'jijpws?; 

Run • the number OF The run . . 

_ 5 -gcORfiELATION distance <SEE TEXT! ^ 

OvnTRACK distance (IN FEET) OF THE 
ISTANCe TO the SOUTHERN/EASTERN Fli 


TMf 

THf 

THg 


^ ^OUNpARIES. 


data 

ELD 


ISISZ • tHE ARRAY CONTAINING THE SCATTEROMETCR DATA 
THE INPUTS FROM UNIT 10 AREt 
IFLD • AN ARRAY CONTAINING THE NUMBERS OF THE SAMPLED 
start ■ Tm|*'oSvnTRACK D^sfANcks 5 f THE FIELDS IN THE LINE 

?!18f : ?SI 

THE OUTPUT variables AREt 

X ■ NUMPES of FEET FRO 

Ic.T : 0 . 

SetCCCCCJJg; 

REAL 


ABOVE 

- distance 

FROM northern field BOUNDARY 

from the mest|rn f|elo boundary 


rFal 


[GZ(S10*|0) fSCA 


S) tYWlDEUS) fSCA 
(SIO) 


(||t|OWSHEAR(lS) *X(20> *Y(20 


2 IFLOdSttlSI 
RT(15)«XWlOE(i 
0RR(5I(»).tSCATY 

READ the input files 

format ( lXt3Il4lXiIl.lXtIl.lX* 

^REAo'(9doltENO»ll>DECORR(M» »SCATX(MI *$CATY(M) . (ISI6Z(M*N) *N»l*10 
FORMAT (9X*F4.0*lXtF6.0*IX*FS.0.6X.10IA) 


100 


».IC3.ILlNE.IRUN*0ECORR(K).SCATXm .SCATY(l) f 
»F4*0.1X.F6.0*1X.FS.O.OX.10U) 


101 

li 


SCATY (H) bSCATY (M) . 630*0 
NUMSfM 


CONTINUE . 
00 2 NbI.IS 


10 


8 

C 


REAO(i0.102.ENO«12l IFLO(NI .START (N) *XWIDE (N) tYWiDE (N) 
F0 RMaT(12*1X.F7.1*!X.F6.1.1X.F6*1) 

NFLOSsN 

M. IS THE COUNTER FOR INPUT FILE 10 
12 M«l 

THIS SECTION LOOKS FOR THE BOUNDARIES IN THE DATA 


DO 3 KbI.NUMB 

if(scatx(k) ,lt.staht(m)igo to 

knOB^TART (Ml .XWIDE (M) 


4 

44 


DO 4 NbI«20 
IF(SC “ 


CATX (U .GT.ENDIGO TO 44 


L*L.1_ 

ARANk*SCATX (U -SCATX (K) 

THIS SECTION SETS UP THE FIELD DIMENSIONS 


NPTSbL-K.I 

Bo.r 


BTWN 

XVbXU 

YVbY>* 

F(IL 

F(!l 

F(IL 


(M) 

_.(M> 
nE.EO.DXWbO.O 
NE.E0.4)XWbO.O 
NC.SE.6)YWbO*0 


8S 

IS 

8S 

6R 

GR 

GR 

GM 

6H 

GH 

GR 

GR 

6R 

GR 

GR 

GR 

Gh 

G'’ 

GR 

GR 

GR 

GR 

^R 

8; 

GR 

GH 

8S 

GR 

GR 

)GR 

88 

GR 

GR 

GR 

83 

GR 

GR 

)GR 

6^ 

GR 

GR 

83 

GR 

GR 

GR 

GR 

GR 

GR 

GR 

GR 

GR 

GR 

GR 

GR 

GR 

GR 

GR 

GR 

GR 

GR 

83 

6h 

GH 

6H 

6H 

6H 

GH 

GH 

GH 

GH 


Figure A-4.- Coinputer listing for program GRID. 


00390 

00400 

00410 

00420 

00430 

00440 

004S0 

00460 

00470 

00480 

00490 

00500 

00510 

00520 

00530 

00540 

00550 

00560 

88li8 

00590 

OU600 

00610 

00620 

00630 

00640 

00650 

00660 

00670 

00680 

00690 

00700 

00710 

00720 

00730 

00740 

00750 

00760 

00770 

00780 

00790 


A-7 


riLCl OAIO 


rORTAAN A 


UKUaiNAL PAGE^ 

OF POOR QUALITY 

conversational monitor system 


i 

c 


IP(ILINE*OT.4)00 to as 

this section is for latituoinal flioht lines 

E06E> I xw IDE CM) -ARANGE ) /2 • 0 
00 5 N»ltNPTS 
MUM^L-N*1 

X(N>«IFIXUBS(XW-(CDOE*NTVNm 
r <N)»!FlxiY«lOE(M»-SCATYCNOH)» 

SMEAR (N)bIFIX(0EC0Rr7numM 

00 A NNaltlO 

6 SCAT(N*NN)sISIGZINUM»NN| 

XTR«L*N 

GTVN«SCATX (L > -SCATX (KtRl 

» 

THIS SECTION IS FOR LONGITUOINAL FLIGHT LINES 


NOM*L-N* 1 

X(N)«IFix(XV«SCATV(NUM) I 
Y (N) •IfIx ( ABSTyw- (E0GE«STWN) 
SMEARlNIsIFlXtOECORHtNUM) ) 

00 A NN«I*10 


)) 


SCAT (Nf NN) «ISIGZ (NUMtNN) 

KTR«L-N 
BTWN«SCATX(L)-SCATX(KTR) 

T CONTINUE 

PREPARE TO WRITE THE OUTPUT FIL' S 
S5 IUNIT«10^N 

WRITE OUT THE FIELOS COMPUTED IN THIS PROGRAM 
00 22 KK«1.NPTS 

22 WRiTE(lUNlT«lU)lCl*IC2*IC3«lFL0(M),ILlNEtIRUNtSMEAR(KK)» 
6 X(KK) fYlKK) t (SCAT (KK*MN)«MN«I,10) 

FORMAT(}Il.U*I2tIX*nfIl*lXtI3*IXflAfIXfl4»lXf.lOlA) 

IUNITbIUNIT*! 

IF(M. 


Ill 


MaM^l 
3 CONTINUE 
33 _ CONTINUE 


EO.NFLOS)GO TO 33 


STO? 

end 




ooaoo 

ooaid 

ooazo 

00030 

OOSAO 

ooaso 


so 


in? 


Figure A-4.- Concluded. 


A-8 


ORluiNAL 

OF POOR QU ■ r 


riLCi ORIO 


cxec 


conversational nonitoh system 





HOUSTON 

CMSLIB fortran 

aoER (PERM 
SK GRID LISTING Al 
INCH (PERM_ 

RHINAL (PERM _ . 

.K il YOATA a (PERM RECFM F L 
,A BNORY* AOATA A (PE«M^RECFM 
DISK L3 ?0AIA A (PERM RECFM F LRECL 


(PERM RECFM FA 
lRE 


PERM RECEM F 


* ' 

ra'tiS-iilST8*-58«‘ ‘ 

>'2rci**ni8t8‘-J8S|‘ * 

°IIoto^2(>Jne * 

50 DISK L12 20ATA 
» 10 LGOTO -(JOnC^ 

^.YPE NOT ENOUGH FILEOEF*! 

i|RR LTYPE ERROR ENCOUNTERED IN DEFINING FILES 

-DONE LTYPE 12 OUTPUT FILES HAVE BEEN DEFINED 
LOAD GRID 



(PERM RECFM F 
(PERM RECFM F 
(PERM RECFM F 
(PERM RECFH F 
(PERM RECFM F 
A (PERM RECFM F 
A (PERM RECFM F 
A (PERM RECFM F 
A (PERM RECFM F 
A (PERM RECFM F 
A (PERM RECFM F 
I IN EXEC 


» 
LRECL 6« 
LRECL G4 
LRECL G4 
LRECL 64 
LRECL 64 
LRECL 64 
LRECL 6a 
LRECL 64 
LRECL 64 
LRECL 64 
LRECL 64 
LRECL 64 


GO 


Figure A-5.- Execute file for program GRID. 
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